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The dynamics of a superconducting quantum point contact biased at subgap voltages is shown to be 
strongly afltected by a microwave electromagnetic field. Interference among a sequence of temporally 
localized, microwave-induced Landau-Zener transitions between current carrying Andreev levels 
results in energy absorption and in an increase of the subgap current by several orders of magnitude. 
The contact is an interferometer in the sense that the current is an oscillatory function of the 
inverse bias voltage. Possible applications to Andreev-level spectroscopy and microwave detection 
are discussed. 
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The classic double-slit interference experiment, where 
two spatially separated trajectories combine to form an 
interference pattern, clearly demonstrates the wavelike 
nature of electron propagation. For a 0-dimensional sys- 
tem, with no spatial structure, a completely analogous 
interference phenomenon may occur between two distinct 
trajectories in the temporal evolution of a quantum sys- 
tem. Such trajectories may appear in the presence of 
temporally localized non-adiabatic perturbations (rather 
than spatially localized slits in a screen) which scatter the 
system from one adiabatically evolving state to another. 
In this Letter we show that this type of interference phe- 
nomenon significantly controls the microscopic dynamics 
of a voltage-biased superconducting quantum point con- 
tact (QPC) subject to microwave irradiation. 

It is well known that the Josephson current in a QPC 
is carried by Andreev bound states localized within the 
contact area. The corresponding energy levels — the 
Andreev levels — lie in the energy gap of the supercon- 
ductor and their positions depend on the change in the 
phase of the condensate across the junction. Hence, the 
Andreev levels will move adiabatically with time within 
the gap if the contact is biased by a voltage V much 
smaller than the gap energy, A. With any (normal) elec- 
tron scattering present in the contact the Andreev levels 
will, however, never cross the Fermi level; instead they 
will oscillate periodically with so that on the average no 
energy is transfered to the QPC and a purely ac current 
will flow through the contact (ac Josephson effect). 

Microwave radiation of large frequency w ~ A repre- 
sents a non-adiabatic perturbation of the QPC system. 
However, if the amplitude of the electromagnetic field 
is sufficiently small, the field will not affect the adiabatic 
dynamics of the system much unless the condition for res- 
onant optical interlevel transitions is fulfilled. Such res- 
onances will only occur at certain moments determined 
by the time evolution of the Andreev level spacing. They 
provide a mechanism for the energy transfer to the sys- 



tem to be nonzero when averaged over time and for a 
finite dc current through the junction. The rate of en- 
ergy transfer is in an essential way determined by the in- 
terference between different scattering events Q], which 
will also lead to oscillatory features in the current-voltage 
characteristics of the QPC. 

Although an example of the more general problem of 
energy level spectroscopy, the spectroscopy of Andreev 
levels has an important specific aspect. Due to their 
ability to carry electric current, detection of optical tran- 
sitions between Andreev levels is possible by means of 
transport measurements. The appearance of a subgap 
current under resonant radiation can furthermore be used 
as a sensitive microwave detector. 
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FIG. 1. Time evolution of Andreev levels (full lines) in the 
energy gap of a gated volt age- biased, single-mode supercon- 
ducting QPC (see inset). A weak microwave field induces 
resonant transitions (wavy lines) between the levels at points 
A and B and the level above the Fermi energy becomes partly 
occupied to an extent determined by interference between the 
two transition amplitudes. Non-adiabatic interactions release 
the energy of quasiparticles in the (partly) occupied Andreev 
level into the continuum at point C, where the Andreev states 
and the continuum merge into each other (represented by 
dashed arrows, see text) and the initial conditions for the 
Andreev level populations are reset (filled and empty circles). 

For an unbiased QPC, the Andreev spectrum of each 
transport mode has the form i?^(0) = ±E{(j)) = ±A[1 — 
-Dsin^(0/2)]^/^, where D is the transparency of the mode 
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and the energy is measured from the Fermi energy . 
With a small bias voltage applied, the levels move along 
the adiabatic trajectories E^{t) = ±E{(j)o + 2eVt/h) in 
energy-time space, as shown in Fig. 1. When the cri- 
terium U^' = ^ 2E^{t)/A for adiabaticity is obeyed, 
the rate of interlevel transitions is exponentially small 
keeping the level populations constant in time |^,^ . The 
presence of a weak electromagnetic field [on the scale of 
E{t)] does not affect the adiabatic level trajectories ex- 
cept for short times close to the resonances at t = tA,B, 
when E{tA,B) = hLo/2. Here the dynamics of the system 
is strongly non-adiabatic with a resonant coupling which 
effectively mixes the adiabatic levels. This is an analog of 
the well known Landau-Zener transition, which describes 
interlevel scattering as a resonance point is passed. In 
our case these transitions give rise to a splitting of the 
quasiparticle trajectory at the points A, B into two paths; 
A1A2B2 and A1B1B2, forming a loop in {E,t) space. 

The resonant scattering opens a channel for energy ab- 
sorption by the system; a populated upper level when 
approaching the edge of the energy gap (at point C in 
Fig. 1) creates real excitations in the continuum spec- 
trum, which carry away the accumulated energy from 
the contact. As a result, the net rate of energy trans- 
fer to the system is finite; it consists of energy absorbed 
both from the electromagnetic field and from the voltage 
source. The confluence of the two adiabatic trajectories 
at B2 (see Fig. 1) gives rise to a strong interference pat- 
tern in the probability for real excitations at the band 
edge, point C. The interference effect is controlled by 
the difference of the phases acquired by the system dur- 
ing propagation along the paths A1A2B2 and AiBiB2- 

For a quantitative discussion we consider a one-mode 
superconducting point contact with arbitrary energy- 
independent transparency D for normal electrons, < 
D < 1. The contact is biased at a small applied volt- 
age eV <ti A, and a high frequency electromagnetic 
field is applied to the gate situated near the contact, 
see inset in Fig. 1. We will describe the evolution of 
the Andreev states with the time-dependent Bogoliubov- 
de Gennes equation for the quasiclassical envelopes 
u±{x,t) of the two-component wave function 'I'(x,t) = 
M+(a;, i)e*'=^^ -t- u_(a;, i)e-*'=^^, 



ihdu/dt = [Ho + (T^Vg{x, t)]u . 



(1) 



In this equation u = (w+, u_) is a four-component vector, 
Hq is the Hamiltonian of the electrons in the electrodes 
of the point contact. 



Hq — —ihvFazTzd/dx + 

A[cos{(f){t)/2)ax + i sin((/)(i)/2)sgna;crj^], 



(2) 



where Ui and denote Pauli matrices in electron-hole 
space and in u± space respectively. The function u, 
which is smooth on the scale of the Fermi wavelength, 
has a discontinuity at the contact which is determined 



by the transfer matrix of the QPC in the normal state 
and is described by the following boundary condition [^j : 
n{+0) = {l/y/D){l-VRTy)u{-Q), R=\-D. 

The gate potential Vg{x,t) = Vi^{x) costot in Eq. (|l|) 
oscillates rapidly in time and the amplitude is assumed 
to be small compared to the Andreev level spacing, 
<C E{t). Under this condition, the system experiences 
an adiabatic evolution at all times except close to the res- 
onances (points A and B in Fig. 1) The duration 6t of 
these resonances is short in the limit eV A compared 
to the period of Josephson oscillations Ty = nh/eV. In- 
deed, the resonant transition occurs if the deviation of the 
interlevel spacing from the resonance value E{t) — hLo/2 = 
E{tA.B)^t does not exceed the quantum mechanical res- 
olution of the energy levels h/5t. From this we can esti- 
mate 5t to 5t ~ [Ti/ EltA^B)]^^"^ Ti/eV. Hence we may 
consider the non-adiabatic dynamics as temporally local- 
ized scattering events. By introducing a linear combina- 
tion u(t) = 6+(t)e'"*/2u+ -I- 6-(t)e-*"*/2u- of the eigen- 
statcs u='= corresponding to the adiabatic Andreev levels 
E^, we can describe the system's evolution through a 
resonance by letting a scattering matrix S connect the 
coefficients Jy^ before and after the splitting points A 
and B. A standard analyses of the Landau-Zener in- 
terlevel transitions (see e.g. [^) gives the scattering ma- 
trix elements at the point A {Sa)++ — {Sa) = 

{Sa)+^ = = where jpp = l-|rp = l-er^ 
is the probability of the Landau-Zener interlevel tran- 
sition. Here 7 = 7r|V^ \'^/{dE/dt) where is the 

matrix element for the interlevel transitions. At the 
splitting point B the scattering matrix reads Sb = S'^. 

The matrix element was calculated for the case of 

a double barrier QPC structure in Ref. For a sin- 
gle barrier junction an analogous calculation gives us, 
V+- = a{L/^o)^/DRVuj sin(0/2), where L is the length 
of the normal region, ^0 is the coherence length, and the 
constant a ~ 1 is determined by the position of the bar- 
rier. We note, that this matrix element is proportional to 
the reflectivity of the junction and therefore the effect un- 
der consideration does not exist in perfectly transparent 
QPC (£» = !), cf. Refs. |,|l0|. 

By introducing the matrix $j ^ = exp (i(Tz$(i, j)). 
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( E{<j^) - ^\ , (3) 



which describes the "ballistic" dynamics of the system 
between the Landau-Zener scattering events, we connect 
the coefficients hr^ at the end of the period of the Joseph- 
son oscillation, (f) — with the coefficients 6^ at the 
beginning of the period, cj) — — 27r, 



^Q,bSb^B,aSa^A,-2tt \ ^_ 



(4) 



The time-averaged current through the junction can be 
directly expressed through these coefficients. 
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The quasiclassical equation for the total time depen- 
dent current at the junction {x = 0) reads I{t) = 
vpiuazTzu), where (..) denotes a scalar product in 4- 
dimensional space. From Eqs. (0J|) it follows that, 



I{t) 



2e 



dt 



dx 



ih 



d{uu) 
dt 



Va 



(5) 



In the static limit, Vg — and 0^0, it equals the usual 
equation / = [2e/K) {dE^ /dip) for the Andreev level cur- 
rent. In the general nonstationary case u is a linear com- 
bination of u='= and we calculate the current averaged 
over the period Ty. Using the normalization condition 
l^+P -I- = 1 and omitting small contributions from 
rapidly oscillating terms, we obtain 
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huj 
~2 



(6) 



The direct current through the contact can be viewed as 
resulting from photon-assisted pair tunneling or equiva- 
lently as being due to the distortion of the ac pair current 
due to the induced interlevel transitions. The magnitude 
of the current is such that the energy absorbed from the 
voltage source, V -Idc, together with the energy absorbed 
from the hf field corresponds to the energy necessary for 
creating a real continuum-state excitation. 

Let us now discuss the boundary condition at (/) = 27rn 
(n is an integer). In the vicinity of these points, the An- 
dreev levels approach the continuum and the adiabatic 
approximation is unsatisfactory, even at small applied 
voltages and weak electromagnetic fields. The duration 
bt of the non-adiabatic interaction between the Andreev 
level and the continuum states can be estimated using the 
same argument as for the microwave-induced Landau- 
Zener scattering. One finds that 6t ~ ^/(Ae^y^)^/^. To 
derive the boundary condition, for example, at point C 
in Fig. 1 , one needs to calculate the transition amplitude 
connecting the states u+(ti) at time t\ tc — St and 
u+{t2) at time > tc + St: (u+(i2) U(t2, ti)u+(ti)). 
Here U(t2,^i) is the exact propagator corresponding to 
the Hamiltonian in Eq. (|l]). It follows from symme- 
try arguments that this amplitude is exactly zero. Both 
the Hamiltonian (||) and the boundary condition for u 
at cc = are invariant under the simultaneous charge- 
and parity inversion described by the unitary operator 
A = PaxTz, where P is the parity operator in x-space. 
This implies that at any time any non-degenerate eigen- 
state of the Hamiltonian is an eigenstate of the symmetry 
operator A with eigenvalue +1 or —1 and that this prop- 
erty persists during the time evolution of the state. In 
particular, the static Andreev state obeys the equation 
Au+((/)) = Au+((/)) at any 0. It follows from Eq. (||) 
that u+ (-(/)) = aa:TyU+{4)) so A(-(/)) = -A{(j)). The im- 
mediate consequence of this property is that the static 
Andreev states correspond to different eigenvalues A at 



(/> < and at > 0, and therefore they are orthog- 
onal ||ll|, even though they belong to the same level. 
Hence the state evolving from the adiabatic state u+(ti) 
is orthogonal to the adiabatic state u+(t2)- As a re- 
sult the probability for an adiabatic Andreev state to be 
"scattered" into a localized state after passing the non- 
adiabatic region is identically zero. In reality, the An- 
dreev state as it approaches the continuum band edge 
decays into the states of the continuum. Such a decay 
corresponds to a delocalization in real space and is the 
mechanism for transferring energy to the reservoir p^ . 

The orthogonality property shown above guarantees 
that the coherent evolution of our system persists dur- 
ing only one period of the Josephson oscillation and that 
the equilibrium population of the Andreev levels is reset 
at each point (j) = 27rn |l3| . This imposes the boundary 
conditions 6+ (27rri+0) = 0, 5^(27rn-|-0) = 1 in the begin- 
ning of each period. Combining this boundary condition 
with Eqs. (0) and (0), we finally get 
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(7) 



where is the phase of the probability amplitude for the 
Landau-Zener transition, which weakly depends on V . 
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FIG. 2. Current vs. inverse voltage from Eq. (7) for a 
biased superconducting QPC irradiated with microwaves of 
frequency uj — 1.52A/R and amplitude corresponding to a 
matrix element |V±| = 0.024A for interlevel transitions. Note 
the cut in the inverse voltage scale. Results of the scat- 
tering approach, Eq. (5) [<>], are close to those obtained by 
numerically solving the BdG-equation (1) with the radiation 
field treated in the resonance approximation (solid line) . The 
close fit means that the scattering picture can be used to re- 
construct the Andreev levels from the period of the current 
oscillations (Andreev level spectroscopy, see text). 

Equation (^) is the basis for presenting the biased QPC 
as a quantum interferometer. There is a clear analogy be- 
tween the QPC interferometer and a standard SQUID in 
that they both rely on the presence of trajectories that 
form a closed loop. In a SQUID, which is used to measure 
magnetic fields, the loop is determined by the device ge- 
ometry; in the QPC the voltage (analog of the magnetic 
field) is well defined while the "geometry" of the loop in 
(E, t) space can be measured. This loop is determined by 
the Andreev-level trajectories in (i?, t) space and is con- 
trolled by the frequency of the external field. This gives 
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us an immediate possibility to reconstruct the phase de- 
pendence of the Andreev levels from the frequency de- 
pendence of the period 11 of oscillations of the current 
versus inverse voltage, see Fig. 2. Indeed, it follows from 
Eq. (|), that 



TTzb 



duj 



(8) 



j=2E 



In order to be able to do interferometry it is necessary 
to keep phase coherence during one period of the Joseph- 
son oscillation. There are three dephasing mechanisms 
that impose limitations in practice: (i) deviations from 
an ideal voltage bias, (ii) microscopic interactions, and 
(iii) radiation induced transitions to continuum states. 
The main source of fluctuations of the applied voltage is 
the ac Josephson effect. According to the RSJ model, 
a fixed voltage across the junction can only be main- 
tained if the ratio between the intrinsic resistance Ri of 
the voltage source and the normal junction resistance Rjsi 
is small. If Ri / Rn ^ 1 the amplitude of the voltage fluc- 
tuations 5V is estimated as 5V ^ {Ri/RM)V . Effects of 
voltage fluctuations on the accumulated phase ^(A, B) 
can be neglected if 5^{A, B) = ^'{A, B)5V < 27r, i.e. if 
eV > {Ri/ Rn)^, which corresponds to a lower limit for 
the bias voltage voltage . 

The dephasing time due to microscopic interactions 
is comparable to the corresponding relaxation time . 
This mechanism of dephasing can be neglected as soon 
as the relaxation time exceeds the Josephson oscillation 
period, 3> h/eV. Taking electron-phonon interaction 
as the leading mechanism of inelastic relaxation, we esti- 
mate Ti to be of the order of the electron-phonon mean 
free time at the critical temperature, Tph(Tc), since the 
large deviations from equilibrium in our case occur in the 
energy interval E < A. This gives jl^ another limitation 
on how small the applied voltage van be, eV > 10~^A. 

The third mechanism of dephasing becomes important 
when the Andreev levels are closer than huj to the con- 
tinuum band edge. One can estimate the correspond- 
ing relaxation time as ^ hA/V/^. For small radia- 
tion amplitudes exceeds Ty, while for optimal am- 
plitudes they are about equal. The effect of the level- 
continuum transitions on the interference oscillations de- 
pends on the frequency. If hui < 2 A /3 the "loop region" 
[E{t) < huj/2] is optically disconnected from the contin- 
uum, and transitions cannot destroy interference. Possi- 
ble level-continuum transitions at times outside the loop 
will only decrease the amplitude of the effect by a factor 
exp{—aV^/eVA), where a < 1 is the relative fraction of 
the period Ty during which transitions to the continuum 
are possible. Accordingly, this factor is of the order of 
unity for the voltages that correspond to the maximum 
amplitude of oscillations. If hui > 2 A/3, the interference 
is impeded by the optical transitions into the continuum, 
and the current oscillations decrease. Still, a nonzero 
average current through the junction will persist. 



The interference effect presented here can also be ap- 
plied for detecting weak electromagnetic signals up to 
the gap frequency. Due to the resonant character of the 
phenomenon, the current response is proportional to the 
ratio between the amplitude of the applied field and the 
applied voltage, I ~ |V±p/Aey. At the same time, for 
common SIS detectors a non-resonant current response 
is proportional to the ratio between the amplitude and 
the frequency of the apphed radiation [0), I ^ |V±/wp, 
i.e. it depends entirely on the parameters of the external 
signal and cannot be improved. 

In conclusion, we have shown that irradiation of a 
voltage biased superconducting QPC at frequencies uj 
A can remove the suppression of subgap dc transport 
through Andreev levels. Due to the resonant nature of 
the photon-induced interlevel scattering the phenomenon 
can be used for sensitive microwave detection. Quantum 
interference among the resonant scattering events can be 
used for microwave spectroscopy of the Andreev levels. 

We would like to thank Dag Winkler for productive dis- 
cussions. Support from the Swedish KVA, SSF, Materials 
consortia 9 and 11, and NFR is greatfuUy acknowledged. 



[1] The role of photoinduced interference effects in a normal 
junction was discussed by L. Gorelik, A. Grincwajg, V. 
Kleiner, R.Shekhter, and M. Jonson, Phys. Rev. Lett. 73, 
2260 (1994); L. Gorelik, F. Maa0, R. Shekhter, and M. 
Jonson, Phys. Rev. Lett. 78, 3169 (1997). 

[2] A. Furusaki and M. Tsukada, Physica B 165-166, 967 
(1990). 

[3] C.W.J.Beenakker, Phys. Rev. Lett. 67, 3836 (1991). 

[4] D. Averin and A. Bardas, Phys. Rev. Lett. 75, 1831 (1995); 

Phys. Rev. B 53, R1705 (1996). 
[5] E.N. Bratus', V.S. Shumeiko, E.V. Bezuglyi, and G. 

Wendin, Phys.Rev. B 55, 12666 (1997). 
[6] Our discussion is also valid for a single dominant mode in 

a multi-mode junction. 
[7] P. G. de Gennes, Superconductivity of Metals and Alloys 

(Addison-Wesley, New York, 1966). 
[8] V.S. Shumeiko, G. Wendin, and E.N. Bratus', Phys. Rev. 

B 48, 13129 (1993); V.S. Shumeiko, E.N. Bratus', and G. 

Wendin, Low Temp. Phys. 23, 181 (1997). 
[9] L.Y. Gorelik, V. S. Shumeiko, R. I. Shekhter, G. Wendin 

and M. Jonson, Phys. Rev. Lett. 75, 1162 (1995). 
[10] U. Gunsenheimer and A. Zaikin, Europhys. Lett. 41, 195 

(1998). 

[11] In fact, Andreev level wave functions are 47r-periodic with 
respect to while the energy of Andreev levels and the 
current are 27r-periodic (cf. explicit equations in Ref. |^). 

[12] This conclusion is in agreement with the analysis of multi- 
ple Andreev reflections in Ref. jsj. 

[13] This result is obvious for a pure SNS contact {D = 1) since 
the two Andreev states involved in the transition at the 
band edge have momenta kp and —kp pointing in opposite 



4 



directions. See also y. 
[14] Shunting by the junction capacitance C is small at relevant 

frequences (~l/5t<A) if C< l/iiivA~ 10~^^F. 
[15] E. Shimshoni and Y. Gefen, Ann. Phys. 210, 16 (1991). 
[16] S. Kaplan, C. Chi, D. Langenberg, J. Chang, S. Jafaney, 

and D. Scalapino, Phys. Rev. B 14, 4854 (1976). 
[17] J.R. Tucker and M. J. Feldman, Rev. Mod. Phys. 57, 1055 

(1985). 



5 



